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Abstract 

The impact of recent measurements of heavy-flavour production in deep inelastic ep scattering and 
in pp collisions on parton distribution functions is studied in a QCD analysis in the fixed-flavour 
number scheme at next-to-leading order. Differential cross sections of charm- and beauty-hadron 
production measured by LHCb are used together with inclusive and heavy-flavour production cross 
sections in deep inelastic scattering at HERA. The heavy-flavour data of the LHCb experiment 
impose additional constraints on the gluon and the sea-quark distributions at low partonic fractions 
X of the proton momentum, down to x ~ 5 x 10~^. This kinematic range is currently not covered 
by other experimental data in perturbative QCD fits. 


1 Introduction 


Understanding the nucleon structure is one of the fundamental tasks of modern particle physics. In 
quantum chromodynamics (QCD), the structure of the nucleon is described by parton distribution 
functions (PDFs), which, in collinear factorisation, represent probability densities to find a parton 
of longitudinal fraction x of the nucleon momentum at a factorisation scale fjf. The scale evolution 
of the PDFs is uniquely predicted by the renormalisation group equations for factorisation |[I1I3- 
The jc-dependence cannot be derived from first principles and must be constrained by experimental 
measurements. The precision of the PDFs is of key importance for interpreting the measurements 
in hadronic collisions. In particular, the uncertainty of the proton PDFs must be significantly 
reduced in order to improve the accuracy of theory predictions for Standard Model (SM) processes 
at the LHC. 

Deep inelastic lepton-proton scattering (DIS) experiments cover a broad range in x and fif. In 
the perturbative regime, a wide x-range of I0~'^ < x ^ 10 ^ is probed by the data of the HI and 
ZEUS experiments at the HERA collider Q. These measurements impose the tightest constraints 
on the existing PDFs. However, additional measurements are necessary for a better flavour separa¬ 
tion and to constrain the kinematic ranges of very small and very high x, where the gluon distribu¬ 
tion is poorly known. A better constraint on the high-x gluon is needed for an accurate description 
of the SM backgrounds in searches for new particle production at high masses or momenta. Signif¬ 
icant reduction of the uncertainty of the low-x gluon distribution is important for studies of parton 
dynamics, non-linear and saturation effects. Furthermore, precision of the gluon distribution at 
low X has implications in physics of atmospheric showers, being crucial for cross-section predic¬ 
tions of high-energy neutrino DIS interaction [jH and for calculations of prompt lepton fluxes in 
the atmosphere [Ell. 



Figure 1: Kinematic range in x for the gluon density covered by measurements at HERA and 
EHCb. Eor the HERA inclusive DIS data, the x range is indicated, where the gluon PDE uncer¬ 
tainties are less than 10% at = 10 GeV^. Eor the EHCb data, the upper (lower) edge of the box 
refers to the indicated upper (lower) end of the rapidity, y, range of the heavy-hadron production. 
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Heavy-flavour measurements of the LHCb Collaboration (bllTl at the LHC probe the very for¬ 
ward range of the heavy-hadron rapidity y and are sensitive to the gluon PDF at low x, as schemat- 
ieally shown in Fig. [B For this illustration, in the ealeulation of the ki nematics of heavy-quark 
production at HERA, the leading order (LO) relation is used for the typical gluon x in boson-gluon 

4,„2 

fusion, X = XBj{l + where xgj denotes the Bjorken scaling variable, mg is the heavy-quark 
mass, and ^ is the virtuality of the exchanged electroweak boson. In the case of heavy-quark 

, J p\+n?Q 

production at LHCb, the LO formula x = e g-, assuming — Q in the parton-parton rest 

frame, is applied. Here, pr and p^ represent the transverse and longitudinal momenta of the heavy 
quark, respectively, and Ep is the proton beam energy. 

Heavy-flavour production in proton-proton collisions at LHC is dominated by the gluon-gluon 
fusion process. Therefore the LHCb measurements of charm @ and beauty [Q production in 
the forward region 2.0 <y < 4.5 probe the gluon distribution at 5 x 10“^ < 10“'^, a region 

which is not accessible with HERA data. Note that the LHCb data are sensitive to the product of 
gluon densities in two non-overlapping low and medium-to-high x ranges, as illustrated in Fig. [T] 
Since the medium range is already well constrained by HERA data, which furthermore bridge the 
gap between the two LHCb ranges, the major impact of the LHCb heavy-flavour measurements 
is expected at 5 x 10^^ ^ 10^"^. The advantage of using heavy-flavour data is that the charm 

and beauty masses provide hard scales for the perturbative QCD expansion all the way down to 
their production threshold. To estimate the impact of the LHCb measurement of charm and beauty 
production on the gluon distribution at low x, these data are included in a QCD analysis together 
with the inclusive DIS (Si and heavy-flavour production [[8l|3 cross sections measured at HERA. 


2 Experimental data used in the QCD analysis 

The main objective of the present QCD analysis is to demonstrate the constraining power of the 
measurements of heavy-flavour production in DIS and pp collisions for the determination of the 
PDFs of the proton. The measurements of charm and beauty production at HERA and LHCb, 
together with the combined HERA inclusive cross-section measurements, are used in a next-to- 
leading order (NLO) perturbative QCD (pQCD) analysis. 

Neutral current (NC) and charged current (CC) inclusive DIS cross sections in ep scattering 
are directly sensitive to the valence and sea-quark distributions and probe the gluon distribution 
through scaling violations [[D. HERA measurements of the NC and CC cross sections in DIS at 
a centre-of-mass energy = 320 GeV have been combined taking into account systematic cor¬ 
relations [I3]|. This combined data set contains the complete information on inclusive DIS cross 
sections published by the HI and ZEUS Collaborations based on data collected in the years 1994- 
2000, and has been used for the determination of the PDF set HERAPDFl.O [|3l. The ki nematic 
range of the NC data is 6 x 10“^ < xbj < 0.65, 0.045 < < 30000 GeV^. The CC cross sec¬ 

tions span the kinematic range of 1.3 x 10"^ < xbj < 0.40 and 300 <Q^< 30000 GeV^. These 
combined NC and CC cross sections represent the basis for all PDF determinations. 

In ep scattering, charm and beauty quarks are produced predominantly in the photon-gluon 
fusion process which provides a direct probe of the gluon distribution in the proton. Measurements 
of open-charm production cross sections in DIS at HERA from the HI and ZEUS Collaborations 
have been combined |[8l . Cross sections for charm production were obtained in the kinematic range 
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of 2.5 < < 2000 GeV^ and 3 x 10 ^ < xbj < 5 x 10“^. The combination method accounts for 

the correlations of the systematic uncertainties among the different data sets. These combined 
measurements were used to improve constraints on the gluon distribution and to determine the 
charm-quark mass (HI. The charm reduced cross sections determined as a function of and xbj 
are used in the present analysis together with all provided details on the systematic correlations. 
In addition, cross sections for the production of b quarks in ep scattering, as measured by the 
ZEUS Collaboration (911 are used in the present analysis. These data correspond to an integrated 
luminosity of 354 pb ^ and cover the kinematic range of 5 < 2^ < 1000 GeV^. The b- and c-quark 
content in the events with at least one jet have been extracted using the invariant mass of charged 
tracks associated with secondary vertices and lifetime information, and the Z7-quark mass has been 
measured (^. In the present analysis, the Z7-quark production data are used mainly to improve 
constraints on the Z^-quark mass. 

For additional constraints on the gluon distribution at low x the differential cross sections of 
charm and beauty production in pp collisions al y/s = 1 TeV from the LHCb experiment are used 
for the first time. The measurement of charm production (6l| is based on data corresponding to an 
integrated luminosity of 15 nb^^. Charm production is identified through the full reconstruction 
of decays of the charmed hadron^ D^, D+, Df and A+. The cross sections are measured 
as a function of the transverse momentum, pj, and rapidity, y, of the reconstructed hadrons. The 
LHCb data on 5-meson production in pp collisions (3 correspond to an integrated luminosity of 
0.36 fb”i. The5+,50 and 5® mesons are reconstructed in exclusive decays mainly involving 7/\|/ 
final states. Correlations between the experimental systematic uncertainties are accounted for as 
described in the original publications. An uncorrelated systematic uncertainty is obtained for each 
distribution by subtracting the correlated uncertainties from the total ones. The 3.5% luminosity 
uncertainty is treated as correlated between the measurements of charm and beauty production. In 
the present analysis, the normalised cross sections, for charm and beauty production are 

calculated from the absolute measurements published by LHCb and are used in the QCD analysis, 
with ^ being the cross section in the center bin, 3 <y < 3.5, of the measured rapidity range 

in each pr bin. The uncorrelated experimental uncertainty on ^ is propagated as a correlated 
uncertainty to the respective complementary rapidity bins. The QCD analysis is performed by 
using both, absolute or normalised, representations of the LHCb measurements, alternatively. 


3 Theoretical predictions for heavy-flavour production 

In the QCD analysis, the experimental measurements are confronted with corresponding theoret¬ 
ical predictions. The theoretical predictions for charm and beauty production in both ep and pp 
collisions are obtained at NLO in the fixed-flavour number scheme (FFNS). This scheme and its 
applicability to HERA measurements is discussed in detail in Ref. (HI and references therein. Pre¬ 
dictions for HERA data are obtained by following the approach of the ABM group at NLO using 
its implementation in OPENQCDRAD ifTOl in the framework of HERAFitter (TTI . The number of 
active flavours is set to Nf = 3, and the renormalisation and factorisation scales (pQCD scales) for 

heavy flavour production are chosen as pr = bf = + 4mg, where mq denotes the pole mass 

'Charge conjugation is always implied for charm and beauty hadrons. 
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of c or quark^ For the light-flavour contributions to the inclusive DIS cross sections, the pQCD 
scales are set io Hr= Hf = Q- 

Theoretical predictions for heavy-quark production mpp collisions are obtained using the mas¬ 
sive NLO calculations [fT2] - [T4]l in the FFNS, also available as part of the Mangano-Nason-Ridolfi 

(MNR) calculations ifTSll . The pQCD scales are chosen as prj = with pQ = ^p'j + mg and 

c b 

A^y being coefficients for c and b quarks, which are discussed in the following. These predictions 
were used successfully for beauty production in pp collisions at the SppS ifT^ and the Teva- 
troij^ [fTTll . They are conceptually very similar to the Frixione-Mangano-Nason-Ridolfi (FMNR) 
predictions lf2TI employed for heavy-flavour photoproduction at HERA If22l . 

The cross-section predictions for heavy-flavoured hadron production not only depend on the 
kinematics of the heavy-flavour production mechanism, but also on the fragmentation of the heavy 
quark into a particular final-state hadron. There is no final-state factorisation scale in the FFNS 
since collinear logarithms of the heavy-quark mass are included in fixed-order perturbation theory. 
The calculations in lfT2l - [T4ll^ describe the production of an on-shell heavy quark. Near the kine¬ 
matic threshold, the transition of the heavy quark into the observed heavy-flavoured hadron can be 
taken into account by multiplying the cross section with the appropriate branching fraction. This 
leads to an excellent description of B- and D-meson production measurements at the Tevatron and 
the LHC from = 0 up to pj ~ 4m [l24ll25]l . The scope of these calculations can be extended 
by convoluting the heavy-quark production cross section with a suitable scale-independent frag¬ 
mentation function describing the hadronisation of the heavy quark. The implementation of the 
convolution is not unique once the quark and hadron masses are taken into account, and leads to 
a potentially pr-dependent modelling uncertainty which is, however, small compared to the scale- 
choice uncertainty at NLO. This fragmentation function is used on a purely phenomenological 
basis, since it does not strictly appear in the context of a factorisation theorem, and therefore it 
has to be extracted from data. It depends on the order of the perturbation series but is generally 
assumed to be otherwise universal. Its main effect is to lower the theoretical predictions at large 
Pt- Typical parametrisations used in the literature are those by Peterson et al. |[ 2 ^ depending on 
one parameter £ and by Kartvelishvili et al. |[27]| depending on one parameter a^- 

For the HERA measurements, the fragmentation functions and their uncertainties are consid¬ 
ered and accounted for in the original publications [[8l|9l . The measurements of LHCb are provided 
as hadron-production cross sections and the fragmentation functions have to be applied explicitly 
in order to use these data in the QCD analysis. In addition, fragmentation fractions describing the 
probability of a quark to fragment into a particular hadron have to be applied. The fragmentation 
fractions for c-flavoured hadrons are taken from [[28ll and for Z^-flavoured hadrons from ll29l . 

So far, no fragmentation measurements were performed in pp collisions. Because of similar¬ 
ities of the c-quark production kinematics at HERA and LHCb, the Kartvelishvili fragmentation 

^The pole mass is used for consistency with the pp predictions, since MS running mass predictions are not available 
for LHCb. 

^Provided that fragmentation and other uncertainties are properly treated IfTSlI . Note that the NLO+NLL (FONLL) 
lfT9l calculations used there, and also used by LHCb ||6l|71, slightly reduce the cross sections at high transverse mo¬ 
menta with respect to the pure NLO calculation, while they are identical at low transverse momenta lfT9l . The Tevatron 
data are well described by FONLL ifTSll . This conclusion is also applicable to the NLO calculations lfT4l used here, 
since these were used as input for the NLO part of the FONLL calculations. The claim in iflTll that the NLO predictions 
undershoot the data while FONLL describes them must thus be attributed to parametrisations beyond the perturbative 
part of the calculations, as illustrated, e.g., in Fig. 7 of ll20l . 
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function GTll with Uk = 4.4 ± 1.7, as obtained from corresponding HERA measurements 
extracted for the NLO FENS scheme, is applied for predictions of the EHCb measurements of 
charm-hadron production. The fragmentation is performed in the laboratory frame by rescaling 
the quark three-momentum with the energy of the produced hadron being calculated using the 
hadron mass. This procedure is used for D+ and D+ mesons, and for baryons. For D®- and 
D^-meson production, the contribution from and D*^ mesons is treated as described in ll^ . 
For beauty production, the value CLk = 11 ± 4 is used for all Z?-flavored hadrons, corresponding to 
measurements at FEP . 

The fragmentation-fraction uncertainties are assigned to the measurements and are treated as 
correlated, while the uncertainties arising from the variations of assumptions on the fragmentation 
functions are treated in the form of variations of the theory predictions in the QCD fit. 


4 Details of the QCD analysis 

The open source QCD fit framework for PDF determination HERAFitter ifTTI . version 1.0.0, is 
used. The partons are evolved by using the QCDNUM program |[34l . The analysis-specific modi¬ 
fications to HERAFitter address the heavy-flavour treatment as follows. The massive fixed-flavour 
number scheme OSlI with the number of flavours Nf = 3 is used for the treatment of heavy-flavour 
contributions. The calculation of one-particle inclusive heavy-quark production cross sections in 
hadron collisions at NFO according to lfT4l is implemented by using original routines from the 
MNR code [l36l. The results agree with those obtained with the original MNR code at a level of 
accuracy below 1%. 

The 3-flavour strong coupling constant in the NFO MS scheme is set to = 0.1059 ± 

0.0005, which corresponds to the world average value of as{mz)^f^^ =0.1185 ±0.0006, using 
two-loop evolution equations Il34ll . 

The range of the inclusive HERA data is restricted to ^ = 3.5 GeV^. The procedure 

for the determination of the PDFs follows the approach used in the HERAPDFl .0 QCD fit Q ■ The 
following independent combinations of parton distributions are chosen in the fit procedure at the 
initial scale of the QCD evolution Qq = 1.4 GeV^: the valence-quark distributions xUy{x), xdy{x), 
the gluon distribution xg(x) and the w-type and J-type anti-quark distributions (which are identical 
to the sea-quark distributions), xU(x), xD{x), where xU(x) = xu{x) and xD(x) = xd{x) +xs{x). At 
the scale Qq, the parton distributions are represented by 


XUy{x) 

= (l-x)*^"'' (l±F„^x^), 

(1) 

xd^{x) 


(2) 

xU(x) 

Cl 

1 

(3) 

xD(x) 


(4) 

xg{x) 

= Ag x^« (1-x)''^-Ag x®« (1-x)*^*. 

(5) 


The normalisation parameters Ad^,Ag are determined by the QCD sum rules, the B parameters 
are responsible for the small-x behaviour of the PDFs, and the parameters C describe the shape of 
the distribution as x —)■ 1. A flexible form for the gluon distribution is adopted with the choice 
of C'g = 25 motivated by the approach of the MSTW group Il37l[^ . The 5-quark distribution is 
defined through x-independent strangeness fraction, fs, of the J-type sea, xs = fsxD at Qq, where 
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fs = O.Sl^QQg as in the analysis of ll38ll . including the recent complementary measurement ll^ . 
Additional constraints and = Ajj(l — fs) are imposed, with xu ^ xd as x ^ 0. The 

analysis is performed by fitting the remaining 13 free parameters in Eqs. ([I]-[5]). 

The PDF parameters are determined in HERAFitter by minimisation of a %^-function tak¬ 
ing into account correlated and uncorrelated uncertainties Il40l of the measurements. Systematic 
uncertainties are assumed to be proportional to the central prediction values, whereas statistical 
uncertainties scale with the square root of the predictions. Correlated uncertainties are treated 
using nuisance parameter representation [|4^ . To minimise biases arising from the likelihood tran¬ 
sition to when the scaling of the errors is applied, a logarithmic correction is added to the 
X^-function iHTl . 

The heavy-quark masses are left free in the fit. They are well constrained by the measurements 
of charm and beauty production in DIS and the fitted values (see Table |2] in the Appendix lAl) are 
consistent with the ones obtained in the corresponding HERA analyses EIHl within the intrinsic 
theoretical systematic uncertainty of the pole mass definition ll42l . 

The QCD analysis is performed twice using either absolute or normalised differential cross 
sections of heavy-flavour production from EHCb measurements, as defined in Section |2l The 
implementation of the theory calculations [|43]| as described in Section |3] allows the pQCD scales, 

c b 

i.e. the parameters , and the values for the pole mass of the heavy quarks to be changed at each 
fit iteration. 

In the QCD analysis using the normalised EHCb measurements, the pQCD scales are fixed 
to Af = Af = I for the central result. The scale dependence is studied by varying the pQCD 
scales independently such that 0.5 < Ay.Af < 2. and A^ are always varied simultaneously. The 
resulting scale dependence is small, since it is largely absorbed by the normalisation, as illustrated 
in Appendix IaI 

In the variant of the fit using the absolute EHCb cross sections, the scale dependence of the 
predicted cross section is the dominant theoretical uncertainty. The same scale choice and variation 
procedure, as applied for the variant of the fit using the normalised EHCb measurements, leads to 
unacceptably high y} values of the respective fits ll43]l . Therefore, the four scales technically are 
treated as independent fully correlated systematic uncertainties for the central result. Since the 
pQCD scales are not physical parameters, the related uncertainties are not obtained from the fit. 
Instead, the effect of the scale choice on the other fitted parameters is evaluated by an independent 
variation of Ay^ in the range 0.5<Ay=A^< 2 with A^ and A(^ as free parameters, or Ay in the 
range 0.25 < A^=A^ < 1 with Ay and A^ being free parameters. For the variation Ay = A^ = 0.5, 
a cut pr > 2 GeV is applied for the charm EHCb data to ensure that the factorisation scale is above 
1 GeV^, since this is technically required in the QCDNUM. This procedure ensures an acceptable 
fit quality for all variations [|43l, as required for a meaningful extraction of the other uncertainties. 
Because of the unconventional scale treatment the fit using absolute cross sections is considered to 
be a cross check. 


5 PDF uncertainties 

The PDF uncertainties are estimated following the approach of HERAPDFl.O [|3l in which exper¬ 
imental, model, and parametrisation uncertainties are taken into account. Experimental uncertain¬ 
ties are evaluated using the Hessian method [|40ll . A tolerance criterion of Ax^ = 1 is adopted for 
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defining the fit uneertainties that originate from the experimental uneertainties of the measurements 
ineluded in the analysis. 

Model uneertainties arise from the variations in the values assumed for imposed on the 
HERA data, whieh is varied in the interval 2.5 < <5.0 GeV^; the fraetion of strange quarks, 

varied in the range 0.23 < fs < 0.50 and the value of the strong eoupling, varied in the range 
0.1054 < < 0.1064. The pQCD seales for heavy-quark produetion in DIS are varied 

simultaneously by a faetor of 2 up and down for both, eharm and beauty. For the fits with the LHCb 
data, the model uneertainties inelude theoretieal uneertainties for the eross seetion predietions for 
heavy-flavoured hadron produetion, arising from variation of the pQCD seales and of the fragmen¬ 
tation parameters, as deseribed in Seetion[3l Uneertainties, arising from these model variations are 
referred to as MNR uneertainties in the following. 

The parametrisation uneertainty is estimated similarly to the HERAPDFl.O proeedure: for 
all parton densities, additional parameters are added one by one in the funetional form of the 
parametrisations in Eqs. (IB©, in a similar way as deseribed in |l3l[8l|9l|. Furthermore, the starting 
seale is varied to = 1.9 GeV^. The parametrisation uneertainty is eonstrueted as an enve¬ 
lope built from the maximal differenees between the PDFs resulting from all the parametrisation 
variations and the eentral fit at eaeh x value. The total PDF uneertainty is obtained by adding 
experimental, model and parametrisation uneertainties in quadrature. 


6 Results 

In Fig.[2l the absolute eross seetions for D^- and R+-meson produetion in pp eollisions are shown 
for one representative rapidity bin and are eompared to the theory predietions as used in the QCD 
analysis. A signifieant seale dependenee is observed. The normalised eross seetions for a repre¬ 
sentative pt bin of the same data set are eompared to the respeetive theory predietions in Fig. [3l 
The advantage of using the normalised eross seetion is a signifieant reduetion of the seale de¬ 
pendenee of the theoretieal predietion, retaining the sensitivity of the eross seetions to the gluon 
distribution. The reduetion of the uneertainty due to seale variation is related to the faet that the 
seale ehoiee affeets mostly the normalisation but only to some extent the shape of heavy-quark 
produetion kinematies, as demonstrated in Fig.O|7]in the Appendixj^ 

The fit quality, represented by the total and partial values of divided by the number of 
degrees of freedom, ndof, for both variants of the QCD analysis is presented in Table [IJ When the 
normalised FHCb measurements are used in the QCD analysis, ^dof is appropriately redueed for 
the respeetive data sets. The fitted parameters are presented in Table |2] in the Appendix lAl 

The resulting gluon, valenee-quark and sea-quark distributions with their total uneertainties are 
presented at = 10 GeV^ in Fig. |4] and eompared to the result of the fit, based on solely HERA 
measurements of inelusive and heavy-flavour DIS. The uneertainties on the gluon and sea-quark 
distributions at low x are signifieantly redueed in both oases, using FHCb absolute or normalised 
heavy-quark produetion eross seetions. In ease of the variant of the fit based on normalised FHCb 
eross seetions, the uneertainties are redueed by more than a faetor of three at x ~ 5 x 10~^, whieh 
is the edge of the sensitivity of the ineluded measurements (Fig. [I]). Consistent results are obtained 
in the fit using the absolute eross seetions, whieh is eonsidered an important eross oheek of the 
self-eonsisteney of the NFO theory deseription. 

The individual eontributions of the experimental, model and parametrisation uneertainties for 
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Figure 2: Data to theory eomparison for a representative subset of the LHCb absolute eross seetions 
for produetion of D° mesons for 3.5<y<4.0 (left) and of 5+ mesons for 3.0 <y< 3.5 (right). 
In the bottom panels the ratios theory/data for the nominal variant of the fit and the seale variations 
are shown. For demonstration purpose, eorrelated shifts for data points obtained in the fit using 
nuisanee parameters are applied to theoretieal predietions. Uneorrelated uneertainties for data 
points are shown as they are resealed in the fit, while total uneertainties are shown as not resealed. 


both oases of using the LHCb measurements are shown in Fig.[5]and oompared to the result of the fit 
using only HERA data. The gluon distribution at low x is oonstrained by the HERA measurements 
mostly via the sum rules and this results in large parametrisation uneertainties. Onoe the LHCb 
measurements are ineluded in the QCD analysis, the gluon distribution is direotly probed and the 
parametrisation dependenoe of the PDF is signifieantly reduoed. 

The main differenoes in the PDF uneertainties between the fits using the absolute and nor¬ 
malised LHCb measurements are oaused by the MNR uneertainties. The variation of the pQCD 
soales in the predietion of the absolute eross seetion of heavy-flavour produetion in pp oollisions 
leads to signifieant ehanges in the normalisation of the eross seetion and represents the dominant 
uneertainty on the PDFs. The variations of the assumption on fragmentation parameters [|43]| result 
in a smaller uneertainty, as eompared to that due to the seale variations. 

In the ease of the PDF fit using the normalised LHCb eross seetions, the MNR uneertainty 
is strongly reduoed, sinoe variations of pQCD soales and of the fragmentation parameters do not 
signifieantly affeot the shape of the y distributions for heavy-flavour produetion. Therefore this is 
eonsidered to be the primary result of this paper, while the oonsisteney between the absolute and 
normalised variants is eonsidered to be an important eross eheok. 
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Figure 3: Data to theory eomparison for a representative subset of the LHCb normalised eross 
seetions for produetion of D® mesons for 2.0 < pr < 3.0 GeV (left) and of 5+ mesons for 3.0 < 
Pt < 3.5 GeV (right). The eentral rapidity bins are fixed to 1 by the definition of the normalised 
eross seetions. In the bottom panels the ratios theory/data for the nominal variant of the fit and the 
seale variations are shown. For demonstration purpose, eorrelated shifts for data points obtained in 
the fit using nuisanee parameters are applied to theoretieal predietions. Uneorrelated uneertainties 
for data points are shown as they are resealed in the fit, while total uneertainties are shown as not 
resealed. 

7 Conclusions 

The sensitivity of heavy-flavour produetion in pp eollisions to the low-x gluon distribution was 
studied in a eomprehensive QCD analysis at NLO. The measurements of c- and h-hadron pro¬ 
duetion eross seetions at the LHCb experiment are ineluded into a PDF fit together with inelusive 
and heavy-flavour produetion measurements in DIS at HERA. Sinee the bulk of the heavy-flavour 
data is elose to the kinematie threshold, the fixed-flavour number seheme at next-to-leading order 
is used for the predietions of heavy-flavour produetion in ep and pp eollisions. A signifieant re- 
duetion of the parametrisation uneertainty of the gluon distribution at very low x is observed, as 
eompared to the result of the PDF fit using only HERA DIS data. 

Two ways of using the LHCb measurements in the fit are studied. Although the absolute differ¬ 
ential eross-seetion measurements eontain more information, the resulting PDFs suffer from large 
theoretieal uneertainty due to unealeulated higher-order eorreetions, estimated by the variation of 
the pQCD seales. By using only the rapidity shape information in the normalised eross seetions 
for the final result, this uneertainty is signifieantly redueed for the PDF extraetion. 

The present analysis has illustrated the high potential of the LHCb measurements to eonstrain 
the gluon distribution at low x, and global PDF fits elearly ean profit from the inelusion of sueh data. 
Preeise measurements of normalised eross seetions of heavy-flavour produetion in the forward 
kinematie range of the LHC therefore have a great potential to further improve the eonstraints on 
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PROSA NLO FFNS fit 


PROSA NLO FFNS fit 



PROSA NLO FFNS fit PROSA NLO FFNS fit 



Figure 4: The gluon (top left), the sea-quark (top right), the w-valence quark (bottom left) and the 
(i-valence quark (bottom right) distributions represented at = 10 GeV^, as obtained in the QCD 
analysis of the HERA only data (light shaded band) and HERA and EHCb measurements and their 
relevant uncertainties. The sea-quark distribution is defined d&'L = 2-{u-{-d-\-s). The results of the 
fit using absolute or normalised EHCb measurements are shown by different hatches. The widths 
of the bands represent the total uncertainties. 
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PROSA NLO FFNS fit: HERA DIS 



PROSA NLO FFNS fit: HERA DIS + LHCb abs. 


PROSA NLO FFNS fit: HERA DIS + LHCb norm. 



Figure 5: Individual contributions to the PDF uncertainties on the gluon distributions, obtained 
in QCD analyses using HERA-only (upper panel), HERA+LHCb absolute (lower panel left) and 
HERA+EHCb normalised (lower panel right) cross sections of heavy-flavour production. 
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Representation of the LHCb measurements 

absolute 

normalised 

Global X^/ndof 

1073 / 1087 

958/994 

Global p-value 

0.61 

0.79 

X^-contribution from correlated uncertainties 

73 

49 

X^-contribution from logarithmic correction 

-129 

48 

Data set 

partial X'^/«dof 

NC DIS HERA I combined e’p 

108/145 

108/145 

NC DIS HERA I combined e+p 

419/379 

419/379 

CC DIS HERA I combined e’p 

26/34 

26/34 

CC DIS HERA I combined e^p 

39/34 

41/34 

CC DIS HERA combined 

78/52 

47/52 

bb DIS ZEUS Vertex 

16/17 

12/17 

LHCb 

68/38 

17/30 

LHCb D+ 

53/37 

18/29 

LHCb D*+ 

50/31 

19/22 

LHCb D+ 

24/28 

11/20 

LHCb A+ 

5.3/6 

4.9/3 

LHCb 5+ 

99/135 

81 /108 

LHCb RO 

66/95 

35/76 

LHCb RO 

78/75 

23/60 


Table 1: The global and partial values for the data sets used in the analysis of HERA and LHCb 
measurements. 


the PDFs. 

In order to fully exploit the additional eonstraints from absolute LHC eharm and beauty cross 
sections, a significant reduction of the theoretical uncertainties, e.g. through threshold resum¬ 
mation and/or (partial) NNLO calculations with codes suitable for a usage in QCD analyses, is 
desirable. 
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A Appendix 


Parameter 

absolute 

normalised 

Bs 

-0.14 ±0.07 

-0.08 ±0.10 

Cg 

6.83±0.31 

5.23 ±0.34 


1.74 ±0.22 

1.29±0.32 


-0.19 ±0.04 

-0.16 ±0.05 

Bu^ 

0.668 ±0.020 

0.649 ±0.021 

Cu, 

4.99 ±0.23 

4.98 ±0.23 

Eu, 

12.2±2.4 

13.5±2.7 

Bd. 

0.93 ±0.09 

0.96 ±0.09 

Cd. 

5.50±0.56 

5.59±0.55 

% 

1.63±0.21 

1.63 ±0.24 


0.173 ±0.007 

0.158 ±0.007 

% 

-0.146 ±0.006 

-0.155 ±0.007 

% 

10.4±2.5 

15.1 ±4.2 

ntc 

1.709 ±0.024 

1.257±0.014 

mb 

4.67 ±0.08 

4.19±0.13 


0.659 ±0.020 

1.0 

Ab 

0.262 ±0.007 

1.0 

K 

0.444 ±0.021 

1.0 

4 

0.335 ±0.024 

1.0 


Table 2: The fitted parameters for the NLO QCD analysis using HERA and LHCb measurements. 
The value of strong eoupling = 0.1059 is used. The listed uneertainties eorrespond to 

those assoeiated to the experimental measurements used in the fit. Uneertainties are not quoted for 
parameters that are fixed. 
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Figure 6: NLO QCD predictions for charm LHCb data with different scale choices for absolute 
(top) and normalised (bottom) cross sections. Lower inlets indicate the ratio of predictions to the 
central scale choice. The predictions are obtained by using the FFNS variant of MSTW 2008 
PDFs Il4^ with Nf = 3; the charm mass is set to = 1.5 GeV. 
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Figure 7: NLO QCD predictions for beauty LHCb data with different scale choices for absolute 
(top) and normalised (bottom) cross sections. Lower inlets indicate the ratio of predictions to 
the central scale choice. The predictions are obtained using the FFNS variant of MSTW 2008 
PDFs Il4^ with Nf = 3; the beauty mass is set to mt, = 4.5 GeV. 
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